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Abstract 

Three-dimensional viscous flow 
analysis is performed for a configuration 
where two crossing and glancing shocks 
interact with a turbulent boundary layer. 
A time marching 3-D full Navier-Stokes 
code, called PARC3D, is used to compute 
the flow field and the solution is 
compared to the experimental data 
obtained at the NASA Lewis Research 
Center 1 ft. X 1 ft. supersonic wind 
tunnel facility. The study is carried out 
as part of the continuing code assessment 
program in support of the Generic 
Hypersonic research at NASA Lewis. 
Detailed comparisons of static pressure 
fields and oil flow patterns is made - 
with the corresponding solution on the 
wall containing the shock/boundary layer 
interaction in an effort to validate the 
code for hypersonic inlet applications. 

Introduction 

The shock wave/boundary-layer 
interaction phenomenon is a common 
occurrence in many of the important 


propulsion components such as supersonic 
inlets, nozzles and supersonic 
combustors. The strong pressure 
gradients, secondary flows and possible 
flow separations associated with these 
interactions can have significant effect 
on the flow through the individual 
components and consequently affect the 
overall performance of the engine. 
Therefore, a thorough understanding of 
these phenomena is essential for an 
efficient design of the propulsion 
components in the supersonic flow regime. 
To simulate these complex flows with a 
reasonable accuracy, a full 3-D Navier- 
Stokes code with a fairly reliable 
turbulence model is required. Therefore, 
it is important that the code selected 
for computation of these flows be 
thoroughly assessed using experimental 
data to verify the code's capability to 
correctly simulate the flow features 
mentioned above. The object of this 
study is to assess the numerical code, 
used to compute such flows, using 
experimental data obtained for a 
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configuration where two crossing and 
glancing shocks interact with a turbulent 
boundary layer. The experiment was 
carried out at the NASA Lewis Research 
Center 1 ft X 1 ft supersonic wind 
tunnel facility. 

The numerical code is a time 
marching full Navier-Stokes code, called 
PARC3D 1 , which solves the full three 
dimensional Reynolds averaged Navier- 
Stokes equations in strong conservation 
form using Beam and Warming approximate 
factorization. The code was originally 
developed as AIR3D by Pulliam and 
Steger 2 , and Pulliam 3 later added the 
Jameson 4 artificial dissipation and 
called the code ARC3D. Cooper 1 adapted 
the code for internal flow in propulsion 
application and named the code PARC3D. 
The PARC3D code uses central differencing 
on a generalized curvilinear coordinate 
system with implicit and explicit second 
and fourth order artificial dissipation. 
To simplify the solution of the block 
pentadiagonal system of discretized 
equations, the block implicit operators 
are diagonalized by decomposing the flux 
Jacobians, resulting in a scalar 
pentadiagonal system. The loss of time 
accuracy from the diagonal ization does 
not affect the spatial accuracy of the 
steady state solution. 3 The turbulence 
model used in the code for this study is 
the Baldwin -Lomax model 5 . 

The PARC3D code has been verified 
previously with experimental data for 3-D 
supersonic and 2-D hypersonic flow 
configurations having flow features 
similar to those of the present 
study 6,7,8 . These studies demonstrated 
the capability of PARC3D to simulate with 


reasonable accuracy the shock/boundary- 
layer interaction phenomenon typically 
present in the supersonic inlet 
configurations. This study is a 
continuation of the code validation 
process being carried out at NASA Lewis 
Research Center for computations of flows 
through supersonic/hypersonic propulsion 
components in support of the Generic 
Hypersonic research program. 

Experimental Configuration 

The crossing shock/boundary-layer 
experiment was configured by using two 
shock generator plates that span the 
tunnel test section. The 20.3 cm long 
shock generators arranged symmetrically 
produce oblique shocks which cross each 
other on the centerline. The interaction 
of these shocks with the naturally 
occurring incoming flow boundary layer on 
the tunnel wall defines the experiment. 
This interaction is shown schematically 
in Fig. 1. The incoming boundary layer 
was observed to be uniform and fully 
turbulent. The thickness of the boundary 
layer varied between 3.0 to 3.5 cm 
depending on the inflow Mach number. 

The experiment was run at different 
shock generator angles of 4, 6, 10 and 12 
degrees for a Mach number varying from 
2.5 to 4.0. However, all the angles were 
not run at each Mach number because the 
model would unstart and produce a normal 
shock ahead of the shock generator for 
some combinations of Mach number and 
angles; for example for Mach 2.5, a shock 
generator angle of 9 deg. is the limiting 
value. The experimental data include 
surface static pressure measurements and 
oil flow visualization on the tunnel 
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wall. Complete details of the 
experimental configuration and 
instrumentation used in the experiment 
can be found in reference 9. 

Computations 

The computations are performed on 
the Numerical Aerodynamic Simulation's 
Cray-2 computer, located at NASA Ames 
Research Center, using a grid size of 151 
X 71 X 71. The computational grid 
showing one each of the streamwise and 
cross planes is shown in Fig. 2. To 
resolve the viscous layers, the grid 
lines are clustered in regions close to 
the walls using hyperbolic tangent 
functions such that the first grid line 
away from the wall is located at a y + of 
approximately 2.0. Since the 

configuration is symmetric in both 
directions in the cross plane, only a 
quarter of the full flow field is 
computed with symmetry imposed in the 
vertical as well as the transverse 
direction along the planes passing 
through the centerline of the wind 
tunnel. The flow field at the inflow 
boundary, which is a plane just ahead of 
the shock generators, is held fixed. 
This flow field consists of an incoming 
2-D turbulent boundary layer, on the 
tunnel wall, which is matched with that 
of the experiment. Since most of the 
flow in this plane including a large part 
of the boundary layer is supersonic, the 
fixed boundary condition at this plane 
does not violate any physics of the flow. 
At the outflow boundary, once again since 
the flow is mostly supersonic, the flow 
variables are extrapolated from inside. 
The outflow boundary is located several 
planes beyond the shock generator 


trailing edges. The distance between the 
trailing edge and the boundary is about 
25% of the length of the shock generator. 

A nonreflective boundary condition, 
using a simple Mach wave extrapolation, 
is applied on that portion of the 
computational boundary extending from the 
trailing edge of the shock generator to 
the outflow boundary plane. This 
boundary condition allows the expansion 
waves from the trailing edge of the shock 
generator to pass through the free 
boundary. A no-slip condition is applied 
on the solid walls. In the present 
computations the turbulence model is 
modified such that when a reversed flow 
region is encontered, the outer layer 
model chooses the maximum of the two eddy 
viscosity values computed by Baldwin- 
Lomax formulation and P. D. Thomas 10 
model. In addition, based on the study of 
Deiwert 11 , the inner layer is replaced 
with the outer model which extends all 
the way upto the wall in this region. In 
the absence of reversed flow regular 
Baldwin-Lomax two-layer formulation is 
used. 

Results and Discussion 

Solutions are obtained for the cases 
corresponding to the experimental runs of 
Mach numbers 3.5 and 4.0. Only two shock 
generator angles, 6 and 10 deg., are 
considered for computation because the 
experimental flow visualization data are 
available for only these two angles for 
a Mach number of 3.5. For the Mach 3.5 
case both 6 and 10 deg. configurations 
are computed, whereas for the Mach 4.0 
case only 6 deg. geometry is analyzed. 
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Figure 3a and 3b show the static 
pressure contours on the tunnel wall for 
the case of 6 deg. shock generator angle 
and Mach numbers of 3.5 and 4.0 compared 
with the experimental data. The 
comparison shows a good agreement between 
the experimental data and the solution. 
The corresponding static pressure 
contours for the 10 deg. shock generator 
angle and Mach 3.5 case are shown in Fig. 
3c in comparison with the experimental 
data. The comparison shows that the 
prediction in general agrees fairly well 
with the experimental data. 

Static pressure distributions along 
the centerline of the tunnel wall for 6 
deg. shock generator angle and Mach 
numbers of 3.5 and 4.0 are presented in 
Fig. 4a and 4b along with the 
corresponding experimental data. The 
comparison shows a very good agreement 
between the solution and the experimental 
data. The centerline pressure 
distribution for the 10 deg. case and 
Mach 3.5 is compared with the 
experimental data in Fig. 4c. Except for 
a slight di agreement upstream of the 
plateau region of the curve, the solution 
agrees fairly well with the experimental 
data. This slight disgreement is due to 
the discrepency in the location and 
extent of the predicted reversed flow 
region compared to that of the 
experiment. The difference between the 
predicted and the experimental reversed 
flow region and its effect on the 
pressure distribution are discussed later 
in this section. 

Figures 5a and 5b show the particle 
traces close to the tunnel wall from the 
prediction and the oil flow pattern from 


the experimental data respectively for 
the 6 deg. and Mach 3.5 configuration. 
The particle traces are restricted to the 
XY-plane to approximately duplicate the 
oil flow pattern on the wall. Since these 
pictures are supposed to show the flow 
features qualitatively, the agreement 
between the solution and the data is 
good. 

The particle traces and oil flow 
pattern corresponding to the 10 deg. case 
are shown in Fig. 6a and 6b. Both the 
solution and experimental data show a 
reversed flow region in the middle of the 
tunnel wall. However, in the prediction, 
the region is located slightly upstream 
and extends over a larger area compared 
to that in the experimental data. In the 
experimental study 9 it was pointed out 
that the presence of oil in the flow 
could change the character of the flow 
and so the extent and the location of the 
reversed flow region could be different 
from that of the flow without the oil. 
Due to the difference in the location and 
extent of the reversed flow region, the 
predicted pressure distribution shows the 
plateau, a characteristic of the reversed 
flow, slightly upstream and larger in 
extent compared to the experimental value 
(see Fig. 4c). The discrepency between 
the prediction and the experimental data 
in the reversed flow region could also be 
partly attributed to the approximations 
in the turbulence model, used in the 
code, which is an algebraic eddy 
viscosity model. To improve the accuracy 
of predictions in the reversed flow 
regions it might be necessary to use 
higher order turbulence models which are 
currently being investigated at NASA 
Lewis Research center under the code 
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validation efforts in the Generic 
Hypersonic reasearch program. 

Summary 

3-D viscous flow solutions have been 
obtained for a configuration where two 
crossing and glancing shocks interact 
with a turbulent boundary layer in order 
to assess the PARC3D code for hypersonic 
inlet applications. The solutions for 
different shock generator angles and Mach 
numbers have been compared to the 
experimental data which include static 
pressure fields and oil flow patterns on 
the wall which contained the interaction. 
The comparison shows that the PARC3D code 
is capable of predicting the complex flow 
features, typically present in supersonic 
and hypersonic inlets, with a reasonably 
good accuracy. To further improve the 
accuracy in computing the reversed flow 
regions, application of higher order 
turbulence models needs to be examined in 
the future code validation studies. 
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Fig. 2 Computational Grid 
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Fig. 4b Centerline Pressure Distibution for 6 deg. Mach 4.0 



Fig. 4c Centerline Pressure Distibution for 10 deg. Mach 3.5 
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Fig. 5a Particle traces (analysis) close to tunnel wall for 6 deg. Mach 3.5 



Fig. 5b Oil flow pattern (experimental) on tunnel wall for 6 deg. Mach 3.5 

11 


QRISWAL PAGE J3 
CF POOIf QUfcUTY 




Fig. 6a Particle traces (analysis) close to tunnel wall for 10 deg. Mach 3.5 



Fig. 6b Oil flow pattern (experimental) on tunnel wall for 10 deg. Mach 3.5 
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